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In the analysis described in our Letter a wrong sign of the cross sec- 
tion correction for the beam energy spread was chosen. Although this error 
only slightly affected the values of the (f) meson mass and the product of its 
branching ratios, it unfortunately resulted in a strong overestimation of the 
meson total width. By the time the error was discovered, various improve- 
ments in the reconstruction and calibration procedures allowed an increase 
in the accuracy of the cross section measurement. The resulting systematic 
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uncertainty of the cross section and branching ratio measurements is 1.7% 
common for all energy scans. 

The revised results of our analysis are presented in Table |l[ Within the 
errors they are consistent with those in the Letter. The new corrected value 
of the meson total width is: 

I\a = 4.280 ± 0.033 ± 0.025 MeV. 



Scan 


a , nb 


m^ MeV/c 2 


B e+e -B K o K o x 10 4 


1 
2 
3 

4 


1365 ±14 ±23 
1342 ± 13 ± 23 
1396 ± 12 ± 24 
1393 ± 12 ± 24 


1019.506 ± 0.030 ± 0.020 
1019.512 ±0.023 ±0.045 
1019.363 ±0.017 ±0.080 
1019.316 ±0.021 ±0.122 


0.967 ±0.010 ±0.016 
0.950 ±0.009 ±0.016 
0.989 ±0.008 ±0.017 
0.986 ±0.008 ±0.017 


Average 


1376 ± 6 ± 23 


1019.483 ±0.011 ±0.025 


0.975 ±0.004 ±0.017 



Table 1: <fi meson parameters obtained in this analysis 
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Abstract 

About 300 000 e + e~ — » (/) — > K^Kg events in the center of mass 
energy range from 984 to 1040 MeV were used for the measurement 
of the 4> meson parameters. The following results have been obtained: 
<r = (1367 ± 15 ± 21) nb, m = (1019.504 ± 0.011 ± 0.033) MeV/c 2 , 
T^ = (4.477 ± 0.036 ± 0.022) MeV, r e+e --B(0 -► K%K%) = (4.364 ± 
0.048 ± 0.065) • 10~ 4 MeV. 
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1 Introduction 

This paper presents a precise determination of the mass, and total and lep- 
tonic widths of the 0, based on one of its dominant decay modes, — > K^K S . 
High precision measurements of the meson parameters provide valuable 
information for various theoretical models describing interactions of light 
quarks. 

The results of this work are based on data collected with the CMD-2 
detector |T], Q] which has been running at the high luminosity collider VEPP- 
2M || since 1992. The data come from four scans of the center of mass 
energy range 2Eb eam from 984 to 1040 MeV. In the first scan, performed in 
1994 with integrated luminosity of 0.18 pb _1 , the resonance depolarization 
method M was used for precise beam energy calibration at each point. Three 
other scans, performed in 1996, corresponding to integrated luminosity of 
2.19 pb _1 , do not have resonance depolarization information. The collected 
integrated luminosity of 2.37 pb _1 corresponds to 3.4xl0 6 uieson decays. 



2 CMD-2 detector 

The CMD-2 detector has been described in detail elsewhere [0], 0. It is a gen- 
eral purpose detector consisting of a drift chamber (DC) with about 250 jjl res- 
olution in the plane transverse to the beam and proportional Z-chamber (ZC) 
used for the trigger, both inside a thin (0.38 X ) superconducting solenoid 
with a field of 1 T. 

The barrel calorimeter placed outside the solenoid consists of 892 Csl 
crystals of 6x6x15 cm 3 size. It covers polar angles from 0.8 to 2.3 radian. 
The energy resolution for the photons in the Csl calorimeter is about 8 % in 
the energy range from 100 to 700 MeV. 

The trigger signal is generated either by the charged trigger based on DC 
and Z-chamber hits || or by the neutral trigger || which takes into account 
the number of clusters detected in the Csl calorimeter as well as the total 
energy deposition. Two independent triggers of CMD-2 can be used to study 
the trigger efficiency 



3 Analysis 

The process — > K®K S was detected using the K s decay into two charged 
pions. Events were selected according to the following conditions: 



• Two oppositely charged tracks were found coming from the vertex clos- 
est to the beam. 

• Both tracks had polar angles 0.95 < #1,2 < vr — 0.95. The distance 
between the beam and the vertex in the R — tp plane is: R ver t < 1.5 
cm. 

• The invariant mass of the two vertex tracks, taken as pions, was 450 < 
Mi nv < 550 MeV/c 2 (see Fig. flja) and the missing momentum satisfied 
the conditions 60 < P mis < (jE^ eam — m^- +40) MeV/c, where m^o = 
497.672 MeV/c 2 is the neutral kaon mass 0. 

• Track momenta were 140 < P±^ < 300 MeV/c and the average mo- 
mentum was 180 < (Pi + -P2V2 < 250 MeV/c as shown in Fig. [L|b. 

At each energy the cross section of K^K S production was calculated 
according to the formula: 



a(e+e- - K° L K°) 



N 



e ■ L ■ (1 + S rad ) ■ (1 + S loss ) ■ B(K° S ^ ^ 



IX 



where N is the number of events; L is the integrated luminosity determined 
from large angle Bhabha events with the help of the procedure described 
in ||; 5 ra d is the radiative correction calculated according to ||; 5i oss is a 
correction for event losses due to decays in flight and nuclear interactions of 
the charged pions; B(K S — > 7r + 7r~) = 0.6861 ± 0.0028 is the branching ratio 
of the decay K s — > 7r + 7r~ from J7|; and e is a product of the reconstruction 
efficiency and geometrical efficiency (acceptance), e = e rec e geom . The accep- 
tance is the probability to detect two pions from the K s decay within a solid 
angle determined by the cuts above. 

To obtain the number of KlK s events, cosmic ray events were removed by 
subtracting a smooth background in the distribution over the Z-coordinate of 
the vertex as demonstrated in Fig. [l|c. A sum of two gaussian functions was 
used as a fitting curve, one gaussian describing K^Kg events and the other 
describing the shape of background. To estimate a systematic uncertainty 
because of the background subtraction, we used an alternative assumption 
that the background was described by a straight line with a slope. The dif- 
ference in the number of events obtained by two methods was 0.3%. The 
distribution over the K s meson decay length after all selections and back- 
ground subtraction is shown in Fig. |l|d. The exponential curve in the Figure 
was calculated using the generally accepted value of the Kg lifetime and 
agrees well with the data. 



10000 



8000 



6000 



4000 



2000 




400 450 500 550 600 150 

M inv ( MeV/c 2 ) 



200 



150 



• 00 



50 



200 250 300 
P.™™, ( MeV/c ) 




2000 



1500 



1000 



500 



^Vm'i i i i i i /'"iffli'-^V-'fiw 

20 -10 10 20 
Z ( cm ) 




1.5 



K s decay length ( cm ) 



Figure 1: Distributions used for the selection of e + e~ — ► <p — ► K^K^ events, 
a - Invariant mass of two tracks; b - Average momentum of two tracks; c - 
Z-coordinate of the vertex; d - Decay length of Kg meson. 
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Figure 2: Distribution of the space 
angle between two tracks. The his- 
togram is simulation, points with er- 
rors are data. 



Figure 3: Distribution of the K® po- 
lar angle. The histogram is simula- 
tion, points with errors are data. 



The background from the — > K + K~ and <fi — > 7r + 7r~7r° decays was 
estimated from Monte Carlo simulation to be less than 0.04%. 

With the above cuts 2.97 x 10 5 K^K S events have been selected. The 
distribution of the space angle between the tracks presented in Fig. |2] shows 
a minimal angle typical for pions from the Kg decay. Figure [3] shows the 
distribution for the K\ polar angle calculated from the momenta of the tracks 
as well as their angles. The shape of the distribution is consistent with that 
expected for the KK. 

About 50 % of the K\ events have a cluster in the Csl calorimeter re- 
sulting from a nuclear interaction. Using the polar and azimuthal angles of 
this cluster as well as the angles of the clusters produced by pions from the 
decay Kg -^ tt + tt~, one can reconstruct the <fi -^ K^Kg event without DC 
information. Such "test" events were used to determine the charged trigger 
efficiency 

To study the reconstruction efficiency, the same events have been used 
with an additional requirement of ZC hits and at least one track in DC to 



reduce background from the 



7r + 7r 7T° decays. 



About 9 000 "test" events were selected. The detailed description of the 
efficiency determination procedures can be found in |10 . 



The acceptance as well as the correction for decays in flight was calculated 
by Monte Carlo simulation. 



Table 1: Efficiencies, corrections and their errors at 2E, 



beam 



1020.0 MeV 



Efficiency 


Value, % 


Stat, error, % 


Syst. error, % 


E-rec 


90.30 


1.51 


0.70 


t-trig 


98.32 


0.93 


0.50 


£geom 


46.87 


0.35 


0.90 


1+Oloss 


95.53 


0.17 


0.10 


1+^rad 


77.67 


0.26 


0.50 



Table 2: Energy, number of events, integrated luminosity and measured cross 
section for the 1994 data. 



N s 


E c . m ., MeV 


Nr l k s 


L, nb- 1 


a, nb 


1 


1010.272 ±0.030 


135± 13 


16.98±0.25 


43.07±2.84 


2 


1017.086 ±0.020 


8487± 95 


80.50±0.80 


626.73±27.92 


3 


1018.136 ±0.018 


4697± 78 


27.68±0.36 


921.22±81.49 


4 


1018.956 ±0.018 


5347± 74 


25.59±0.34 


1327.92±58.83 


5 


1019.214 ±0.020 


7610± 89 


27.23±1.23 


1325.12±92.42 


6 


1019.986 ±0.020 


9372±103 


35.44±0.39 


1341.46±52.03 


7 


1020.128 ±0.020 


3725± 62 


16.00±0.20 


1164.35±91.97 


8 


1021.850 ±0.022 


1900± 46 


17.34±0.27 


615.33±55.78 


9 


1023.972 ± 0.020 


2645± 90 


48.16±0.45 


277.42±48.23 



Typical values of the efficiencies and corrections are presented in Table p] 
for 2E beam = 1020.0 MeV (the meson peak). 

For the 1994 data the beam energy at each point was measured by the 
resonance depolarization method Q. For the 1996 data the beam energy 
was determined from the momenta of charged kaons and the analysis of the 
magnetic field of the collider M. 

Tables ^,0 show the number of selected events, integrated luminosity and 
measured cross section at each energy. 



4 Data fits and <f> meson parameters 

The experimental points were fit with a function which includes the contri- 
butions of the p, u, <p as well as higher resonances u;(1420), p(1450), u;(1600), 
0(1680), p(1700) (below referred to as u/, p', u>", 4>', p" respectively) and the 



Table 3: Energy, number of events, integrated luminosity and measured cross 
section for the 1996 data. 



No 


E c . m ., MeV 


Nr l k s 


L, nb- 1 


a, nb 






Scan 


-2 




1 


1004.25±0.36 


77± 37 


31.23±0.30 


11.83± 4.26 


2 


1010.53±0.11 


216± 19 


21.62±0.27 


49.24±3.84 


3 


1016.36±0.04 


2265± 50 


29.13±0.33 


398.06±12.26 


4 


1017.19±0.03 


6552± 86 


55.72±0.40 


603.50±12.11 


5 


1018.06±0.03 


8214± 94 


47.04±0.36 


903.80±22.06 


6 


1018.99±0.02 


15755±127 


63.79±0.42 


1295.37±22.60 


7 


1020.00±0.02 


18403±137 


71.46±0.45 


1228.57±18.92 


8 


1020.96±0.02 


10855±108 


52.13±0.38 


926.26±19.72 


9 


1021.88±0.05 


2800± 55 


19.71±0.24 


617.27±26.36 


10 


1027.69±0.09 


862± 70 


21.33±0.29 


136.53±14.05 


11 


1033.63±0.36 


492± 84 


28.76±0.29 


52.61±14.17 


12 


1039.48±0.36 


285± 19 


22.89±0.26 


34.86±4.98 






Scan 


-3 




1 


1004.64±0.36 


39± 28 


42.60±0.45 


4.55± 2.55 


2 


1011.30±0.06 


640± 45 


55.15±0.90 


62.60± 3.62 


3 


1015.99±0.03 


4570± 86 


61.59±0.41 


378.37± 7.97 


4 


1016.93±0.02 


7541±106 


67.59±0.43 


573.30±10.29 


5 


1017.91±0.02 


17836±164 


100.02±0.52 


924.53±11.52 


6 


1019.04±0.02 


27070±190 


103.52±0.53 


1345.42±12.95 


7 


1019.95±0.02 


25681±183 


98.95±0.52 


1296.68±14.12 


8 


1020.86±0.02 


18433±158 


91.79±0.51 


950.55±12.52 


9 


1021.74±0.03 


8685±110 


58.36±0.41 


664.54±12.07 


10 


1022.67±0.05 


3326± 69 


39.55±0.35 


458.51±20.41 


11 


1028.36±0.09 


845± 44 


25.13±0.27 


110.09± 8.24 


12 


1034.06±0.36 


487± 41 


28.44±0.28 


50.00± 8.38 






Scan 


-4 




1 


1004.19±0.36 


46± 10 


25.08±0.26 


11.25± 2.36 


2 


1011.30±0.07 


382± 37 


32.12±0.30 


66.11± 9.68 


3 


1015.91±0.03 


2478± 61 


33.16±0.30 


391.00±18.90 


4 


1016.94±0.02 


6735± 95 


61.84±0.41 


600.76±18.22 


5 


1017.92±0.02 


14703±140 


79.37±0.47 


989.72±18.87 


6 


1018.76±0.02 


21322±163 


84.61±0.48 


1287.46±19.00 


7 


1019.68±0.02 


22823±167 


86.89±0.49 


1356.42±20.34 


8 


1020.68±0.03 


8800±105 


40.75±0.34 


1058.66±29.39 


9 


1021.60±0.03 


6237± 91 


41.57±0.34 


699.48±24.72 


10 


1022.59±0.04 


4072± 74 


38.15±0.32 


465.68±19.25 


11 


1028.41±0.07 


961± 45 


30.77±0.31 

1 H 


104.97±12.18 



nonresonant background: 

v K (s) = ^ 2J ^jft ■ \A bkg + A p + A u + A^ + Ay + A p , + A*,, + A# + A p „f 

where s is the center of mass energy squared, &2K°(s) = (s/4 — m 2 K0 ) z ^ 2 is 
a phase space factor for two kaons, A^g is a constant amplitude of nonres- 
onant background. Ay = {9v^9vkk) ' Ay are vector meson amplitudes and 
Ay = (s — my + iy/sT v (s))~ 1 are vector meson propagators with the energy- 



dependence of the width as in JIT]: 



TJs) = r 



mfaAs) 



r ( \ r (n m l®M R gV^yW y^a^s) 

J- oA-SJ — 1 a; -D27T - r } 5T + - t5 7r°7^ -, 7TT T -D37T T / ^ 

V s$ 2 ,K) $aK) m w *3x(mi 

1 <A S / _ l <t> \ & K+ K- —^ 7 oT + a K h K s — T / 5T + 

+B ®r n (s) | B yfi$ 3 M 
,n $ rn (ml) 37T m^ 3n (m 2 ^ 

where Bf are the branching ratios of the major decay modes of u> and (f), and 
$/(s) are corresponding phase space factors. For the radiative decay into P7 
where P = 77(71"°), $p 7 (s) = (\/s(l — mp/s)/2) 3 . For the $37r(s) calculation 



the model assuming the decay V — > p7r — > 37r was used |12|]. The masses and 
widths of the higher resonances were taken from [[7). 

The constants 9v^9vkk were calculated using the experimental values of 
the electronic widths [|7| and assuming SU(3) relations with ideal mixing. 
The coupling constants of the higher resonances were parameters of the fit. 
Relative phases between p and p'(p"), oj and ui'{ui") as well as between and 
(f)' were taken equal to tc according to |l3| . 

The cross section can be rewritten in terms of do - the cross section in 
the resonance peak for the — > K^K^ decay mode: 

<tk{b) = (To- , 2 s -— hf \ - Ap + -Au + Acj) + 

®2K{m^) s b ' 2 3 

9p'l9p'KK a ; 9p"-f9p"KK a ;; 

9ch9<pKK 9(f>~i9<pKK 

9u'l9u'KK a 1 _ 9u"~/9w"KK a // 

9<h9<pKK 9<l>i9<f>KK 

_9<t>'l94,'KK A J/ _|_ ^6^9 1 2 



9<j>l9<j>KK 9<t>l94>KK 
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Figure 4: Experimental cross section 
and meson excitation curve in the 
channel e+e~ -► <p -► ^^ (all 
data ). 



Figure 5: Residuals between the 
measured and fitted cross sections 
versus center of mass energy in each 
scan 



The non-resonant background amplitude can be written as: 



.4 



bkg 



IO~bkg 



1 



(T m (t> T q 



9<t>~i9<t>KKi 



where a^g is the cross section of the non-resonant background. 

The values of a^g as well as coupling constants of the higher resonances 
obtained from the fit are consistent with zero and were fixed at zero to 
determine the meson parameters. The determination of the beam energy 
H showed that the energy scale of VEPP-2M had a slope of about 



dE 



scale 



dE-, 



0.005, 



beam 



where E sca i e is the beam energy determined and E\, eam is the real beam energy 
of the collider measured by the depolarization method. This slope results in 
a systematic uncertainty of the <p meson width » F scn ' e -Ys ~ 0.022 MeV which 
is smaller than the statistical error of the width value. It allows to consider 
the width as a common parameter for all data fits, while the masses and peak 
cross sections were different for each scan. As it was found from the fit: 



4.477 ±0.036 ±0.022 MeV, 



Table 4: <ft meson parameters obtained in the experiment. 



Scan 


0o, nb 


nty, MeV/c 2 


r e+e -Bxl0 4 , MeV 


1 
2 
3 

4 


1370 ± 26 ± 33 
1332 ±12 ±32 
1374 ± 9 ± 33 
1396 ± 12 ± 34 


1019.487 ±0.064 ±0.020 
1019.558 ±0.027 ±0.045 
1019.444 ±0.014 ±0.080 
1019.399 ±0.027 ±0.122 


4.374 ±0.086 ±0.105 
4.253 ±0.041 ±0.102 
4.386 ±0.027 ±0.105 
4.454 ±0.043 ±0.107 


Average 


1367 ±15 ±21 


1019.504 ±0.011 ±0.033 


4.364 ±0.048 ±0.065 



where the first error is statistical and the second is systematic. The values of 
the meson mass and peak cross section obtained from the fit are presented in 
Table |j. The experimental data together with the fitting curve are presented 
in Fig. [| where the experimental points of individual scans were shifted in 
accordance with the difference between the average values of the mass and 
cross section in the peak and their values in individual scans. The residuals 
between the measured and fitted cross section versus center of mass energy 
are shown in Fig. |5| for each scan. 

The systematic error of the mass from scan 1 in which the beam energy 
was measured by the resonance depolarization method is determined by the 
precision of this technique. Systematic errors for the mass for scans 2-4 
were estimated from the difference between energy values obtained by two 
methods of the beam energy determination mentioned above. For the sec- 
ond scan it was 0.045 MeV and increased to 0.080 MeV and 0.122 MeV for 
scans 3 and 4 respectively. The systematic error of <To consists of two parts. 
The first one is common for all four scans and comes from the uncertainty 
of the radiative corrections to the cross section of Bhabha events used for 
the luminosity determination (~ 1.5%). The second part is about 1.9% and 
comes from various factors: 



K° T K°c 



background subtraction 

reconstruction efficiency 

trigger efficiency 

radiative corrections for the e + e _ — > A^iig process 

correction for decays in flight and nuclear interactions 

solid angle uncertainty 

uncertainty in beam energy spread 

selection criteria and efficiencies for Bhabha events 



- 0.3% 

- 0.7% 

- 0.5% 

- 0.5% 

- 0.1% 

- 0.9% 

- 0.2% 

- 1.3% 



This part of the systematic error of <7o is numerically the same (1.9%) for 
all scans, although it was independently estimated for each scan. Results for 



10 



Table 5: Results of ctq measurements by various groups. 



Group 


a {e + e-^<f>^K° L K s ), nb 


Reference 


OLYA,1978 

OLYA,1984 

ND, 1991 

SND, 1997 


1400 ± 90 
1440 ± 60 
1430 ±143 
1360 ± 46 


m 

[18] 


CMD-2, 1994,1996 


1367 ± 26 


This work 



the (f> meson mass from different scans are also independent. Therefore we can 
average results from different scans both for the cross section in the peak and 
mass. To average data, we used a standard weighted least-squares procedure. 
For calculations of the weights statistical and systematic errors were added in 
quadrature. In case of <jq only the variable part of the systematic error equal 
to 1.9% was added and the resulting error was combined with the common 
error of 1.5%. The average values of the meson parameters are shown in 
the last line of Table f|. 

The value of the (f) meson mass obtained in this work is the most precise 
measurement of this parameter among e + e~ experiments and is consistent 
with all previous results. There is only one fixed target experiment which 
gave a more precise value of the <fi mass |TJj], differing from ours by 2.6 
standard deviations. 

The value of the meson width obtained in our work is consistent with 
all previous results and is the most precise. Note that its precision surpasses 
the world average quoted by [0 , which is based on the results of nine previous 
measurements. 

The measurement of o"o(0 — > KgK^) has been performed by a number of 
groups with the results presented in Table [|. 

One can see that our value o"o(0 — > K^Kg) = (1367 ± 26) nb does not 
contradict these measurements and is the most precise. 

The cross section in the peak obtained in our experiment is related to 
the product T e + e --B((f> — ► K°K S ). To obtain this value, special fits with this 
product as a free parameter have been performed separately for each scan 
(see the last column of Table |j). The averaging procedure similar to that for 
ctq gives the following result: 



B( 



KlKg) 



(4.364 ± 0.048 ± 0.065) • 10~ 4 MeV. 



The obtained value of T e + e - ■ B((j) —> K^Kg) is the most precise direct 
measurement of this product. 
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Using r e + e - from other experiments, one can obtain B(<ft —>■ K^Kg). 
For example, for T e+e - = (1.32 ± 0.04) keV from 0, 5(0 -► K° L K° S ) = 
0.329 ± 0.006 ± 0.010 has been obtained. 

Alternatively, taking B{(f> —> K^K^) from other works, T e + e - can be 
calculated. For B(<p -> K°K°) = 0.331 ± 0.009 from 0, we obtain for the 
leptonic width Y^^, ee = (1.32 ± 0.02 ± 0.04) keV. The first number is our 
experimental error (both statistical and systematic), while the second one 
describes the uncertainty of the values from H. 



5 Conclusion 

Using the CMD-2 data sample of 2.97xl0 5 — > K^Kg events with recon- 
structed K s — > 7t + tt~ decay, the following values of the meson parameters 
have been obtained: 

a = (1367 ±15 ±21) nb, 

m^ = (1019.504 ±0.011 ±0.033) Me V/c 2 , 

r = (4.477 ± 0.036 ± 0.022) MeV, 

r e+e - • B(<p -> K° L K° S ) = (4.364 ± 0.048 ± 0.065) • 10" 4 MeV. 

These results are more precise than the corresponding measurements from 
any other e + e~ experiment, and the value of the width is more precise than 
all previous measurements. 
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